Bipartite entanglement entropy in fractional quantum Hall states 
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We present a detailed analysis of bipartite entanglement entropies in fractional quantum Hall 
(FQH) states, considering both abelian (Laughlin) and non-abelian (Moore-Read) states. We derive 
upper bounds for the entanglement between two subsets of the particles making up the state. We 
also consider the entanglement between spatial regions supporting a FQH state. Using the latter, 
we show how the so-called topological entanglement entropy 7 of a FQH state can be extracted from 
wavefunctions for a limited number of particles. 
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I. INTRODUCTION 

The fractional quantum Hall (FQH) states have long 
fascinated the condensed-matter community due to their 
remarkable transport properties and the exotic nature 
of their quasiparticle excitations. It is in the con- 
text of FQH states that the notion of topological or- 
der in gapped two-dimensional states first arose.™ Re- 
cently there has been enhanced interest in FQH states 
with non-abelian statistics^"^^ due to the possibility of 
implementing quantum computation schemes topologi- 
cally protected from decoherencei^ The unusual features 
of FQH states have been notoriously difficult to charac- 
terize using traditional condensed-matter concepts such 
as local order parameters and n-point correlation func- 
tions. 

In a separate development, recent years have seen 
growing understanding that entanglement measures bor- 
rowed from the discipline of quantum information can 
be useful in probing global features of quantum many- 
particle statesi^iii^ii^ It is thus natural to ask what fea- 
tures of FQH states can be characterized by entangle- 
ment measures. In a recent short report^^ three of the 
present authors have shown that one such entanglement 
measure, the bipartite entanglement entropy, indeed elu- 
cidates the subtle correlations and topological order in 
the simplest FQH states, those in the so-called Laughlin 
sequence. The bipartite entanglement entropy is defined 
by partitioning the system under question into two blocks 
A and B, and using the reduced density matrix of one 
part (e.g., pA — t^^B P obtained by tracing over B de- 
grees of freedom) to calculate the von Neumann entropy 

Sa = - trj£A In pa] ■ 

In Ref. numerical calculations of the entanglement 
entropy between two spatial regions allowed us to extract 
from the Laughlin wavefunctions the so-called topological 
entanglement entropy, a concept introduced in Ref. llllll2l . 
[For brevity we write topological entropy where no con- 
fusion can arise.] In addition, we provided results on 
the entanglement entropy between subsets of the parti- 
cles making up the state. We showed that such particle 



entanglements are bounded by expressions that manifest 
the exclusion statistics in the Laughlin states. 

In this article, we present a systematic discussion of 
bipartite entanglement entropies for FQH states, elab- 
orating on our results in Ref. [Tol . In addition to the 
abelian Laughlin (L) states, we consider a series of non- 
abelian FQH states: the Moore-Read (MR) (or pfafRan) 
states. ^•-'^•^'^'^ In planar geometry, the respective wavefunc- 
tions are given by 



*mr({^.}) - Pf (^—) Y[i^^ - ^jO^e-^- 1^* 



with Pf denoting the antisymmetric PfafRan symbol. We 
shall here consider these same states in spherical geome- 
try, so as to eliminate boundary effects. 

For both series of states, we derive upper bounds 
^bound pai-^jcle entanglement entropies. A marked 
difference between the m — 3 Laughlin state and the 
m — 2 Moore-Read states is that in the latter the lead- 
ing correlations have a 3-body nature, whereas those in 
the Laughlin states are 2-body effects. This difference is 
nicely manifested in the leading terms of a 1/N expan- 
sion of the upper bounds S*^"""*^, which are given by (see 
section HI below for details) 

m — 3 Laughlin state : 
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Another marked difference between the abelian and 
non-abelian states is in the value for the topological en- 
tropy 7. Comparing the Laughlin and Moore- Read states 
at the same filling fraction v — 1/m we have 



7mr 



7l = In y/m , 
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the difference being due to the non-abehan nature of the 
Moore- Read states (see section IV for details) . In this ar- 
ticle we extract values for 7 directly from wavefunctions 
for a limited number of particles in spherical geometry 
(up to = 10 for the m — 3 Laughlin state and up to 
= 18 for the m — 2 Moore-Read state), finding val- 
ues that are consistent with the expected result. These 
results illustrate how the entanglement entropy can be 
used in diagnosing the topological order for a FQH state 
that is only known in the form of wavefunctions for a lim- 
ited number of particles, as is often the case in numerical 
studies. 

In Section |TT] we give a general discussion of possible 
bipartite entanglement measures in itinerant many body 
systems, carefully distinguishing between particle and 
spatial partitioning schemes. In Section IIIII we present 
analytical and numerical results for particle entangle- 
ment in the Laughlin and Moore-Read FQH states. In 
this, the eigenvalue distribution of the reduced density 
matrix /5„^ plays a central role. In subsection IIII Fl we 
relate the eigenvalue distribution for PnA=2 to the two- 
particle correlation function 172 ■ In Section IIVI we dis- 
cuss spatial partitioning, paying particular attention to 
the numerical procedure followed in extracting the topo- 
logical entropy 7. 

While most results in this article are for fermionic FQH 
states (meaning m odd in the Laughlin sequence and m 
even for the Moore-Read states), we briefly comment on 
bosonic states in subsection IIII CI 



II. PARTITION CHOICES ON THE SPHERE 

The entanglement entropy, being a bipartite measure 
of entanglement, depends on the particular partitions be- 
ing considered. Obviously, a many-particle system can be 
partitioned in many ways. Rather than asking which par- 
tition is the "correct" one, we find it more useful to ask 
what information one can extract from various kinds of 
partitioning. Accordingly, we have partitioned both the 
spatial degrees of freedom, and the particles themselves. 
We find that both schemes are useful, for revealing dis- 
tinct features of the many-particle state. 

The kinds of partitioning one is able to study depend 
on the available degrees of freedom. Our calculations are 
all performed for FQH states in a spherical geometryi^i^ 
In this representation the fermions are placed on a sphere 
containing a magnetic monopole. The magnetic orbitals 
of the relevant Landau level are then represented as an- 
gular momentum orbitals; the total angular momentum 
is half the number of flux quanta, L = ^N^. The A^^ -I- 1 
orbitals are labeled either / = to A'^ or Lz = ~L to +L. 
For A^ particles at fractional filling v — 1/m, one finds 
the interesting FQH states for A^^ — mN — S, where S is 
a finite-size shift. The Laughlin states appear at S* = m 
while for the Moore-Read states S = m+1. The "filling" 
acquires the usual meaning 1/ = N /N^ only in the ther- 
modynamic limit. The orbitals are each localized around 



a "circle of latitude" on the sphere, with the / = orbital 
localized near one "pole." 

Since the FQH wavefunctions on a sphere are obtained 
in terms of orbital occupancies, one can either partition 
orbitals or partition particles. Because of the spatial ar- 
rangement of the orbitals, partitioning orbitals is in fact 
equivalent to partitioning spatial regions. The difference 
between spatial and particle partitioning has not been 
stressed in the literature because the most common sys- 
tems studied (in the context of entanglement entropies 
in many-particle states) are spin models, for which there 
is no such distinction. In the cases of itinerant particles 
where there is a difference, the common default scheme 
has been spatial partitioning. In particular, conformal 
field theory results on entanglement scaling^i^^ and the 
distinction between gapless and gapped states observed 
in entanglement scaling^ actually pertain to the block- 
ing of space rather than the particles or spins them- 
selves. The definition of the topological entropy for two- 
dimensional topologically ordered states is also based on 
the entanglement entropy between spatial blocksiiiii^ 

In previous work,^" three of the present authors stud- 
ied the entanglement entropy between subsets of particles 
making up a Laughlin FQH state. We presented upper 
bounds and gave an interpretation in terms of exclusion 
statistics. In this paper we extend these results to the 
Moore-Read states, where the exclusion effects are more 
intricate. We refer to Refs. [iGlfTTlflsl for other studies of 
particle entanglement properties. 

For orbital or spatial partitioning, we define block A to 
be the first I a orbitals, extending spatially from one pole 
of the sphere out to some latitude. In the thermodynamic 
limit, this is equivalent to choosing a disk-shaped block 
A within an infinite planar system. In this limit, since 
each orbital I is associated with a wavefunction of the 
form usual complex coordinate language, a 

disk with I a orbitals has radius cx \/Ta. 

The spatial arrangement of the orbitals constrains us 
to either disc or ring-shaped spatial regions as the A par- 
tition. Since the orbital indexes do not give us access to 
the full two-dimensional degrees of freedom, it is not pos- 
sible to experiment with the various kinds of topologically 
nontrivial partitions suggested by Preskill and Kitaevj^ 
Levin and Wen^i^ and Furukawa and Misguichji^ How- 
ever, as we have reported previously, the spherical ge- 
ometry is sufficient to probe the topological entropy of 
FQH wavefunctions. 

Comparing particle and spatial entanglement en- 
tropies, we remark that the effect of correlations is oppo- 
site between the two cases. For particle partitioning, the 
maximal entropy is realized for uncorrelated fermions; 
correlations tend to lower Sa from the fermion bound S^. 
Spatial entanglement, on the other hand, is entirely due 
to correlations. The point is illustrated by considering 
the m — \ Laughlin state, where the fermions are uncor- 
related. Here the particle entanglement entropy equals 
S"^, while the spatial entanglement entropy vanishes. 
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III. ENTANGLEMENT FOR PARTICLE 
PARTITIONING 

In this section we provide close upper bounds to the 
entropy of entanglement between ua particles of the state 
and the remaining ns — N ~ ua particles. We also dis- 
cuss the n^-particle reduced density matrices that 
arise in this context. 



A. Multiplet structure and particle entropy bounds 

For FQH states on a sphere, the n^i-particle reduced 
density matrices commute with the total angular 
momentum operators L^^ and L^^ of the selected ua 
particles. This implies that the eigenvalues of are 
organized in a multiplet structure of the corresponding 
SU{2) algebra: an eigenvalue for total angular momen- 
tum will be (2L„^ -I- l)-fold degenerate. 

For HA = 2 fermions, each having angular momentum 
L, the 2-particle states have total angular momenta L2 — 
2L—1, 2L— 3, . . ., 1 (0), for L integer (half-integer), giving 
a total number of {2L + l)(2_L)/2 states. A naive upper 
bound to the entanglement entropy is thus 



< ln[{2L + l){2L)/2] . 



(1) 



Inspecting the explicit structure of the fermionic Laugh- 
lin states with m = 3, 5, . . ., one finds that the eigenval- 
ues corresponding to 2-particle states with L2 = 2L — 1, 
2L — 3, . . ., 2L — (m — 2) all vanish. The reason is that 
the correlations in the Laughlin states are such that par- 
ticles cannot come too close together. For example, if 
a first fermion occupies the I = orbital, localized near 
the north pole, the Laughlin wavefunction has zero am- 
plitude for finding a second fermion in orbitals 1 = 1, 
I = 2, . . ., I = m — 1. The highest possible value of the 
angular momentum of the two fermions combined is thus 
L2 = L + (L — to). The remaining number of non-zero 
eigenvalues is (2L -I- (2 — m)){2L + (1 — m))/2, leading to 
an improved bound on the entropy S'„^=2 

SnA=2 < In [(2i + (2 - m)){2L + (1 - to))/2] (2) 

with 2L = m{N — 1) as before. For ha > 2, the multiplet 
structures are more complicated and we need to resort to 
a different method for finding a non-trivial upper bound 
to the particle entropy. In the next subsection we give a 
general derivation for both the Laughlin and the Moore- 
Read series of fermionic FQH states. 



B. Upper bounds for fermionic states 

For N fermionic particles, ua particles in the A block, 
and the total number of orbitals given by iV^-l- 1 = 2L+1, 
fermionic statistics lead to an obvious upper limit Sa to 
the entropy Sa 



Sa<SI= In 



1 



UA 



(3) 



In the FQH states the correlations are such that the par- 
ticles avoid each other and the entropy is further reduced. 
To obtain a handle on this, one may reason as follows. 
The model FQH states in the Laughlin and Moore-Read 
series can be characterized as zero-energy eigenstates of 
a Hamiltonian penalizing pairs and/or triplets of parti- 
cles coming to the same position. After tracing out the 
coordinates for the B set, the dependence on those in the 
A set is such that one still has a zero-energy eigenstate. 
However, the number of orbitals available to the A parti- 
cles is larger than what is needed to make the model FQH 
state in the A sector, and one instead has a certain num- 
ber of quasi-holes on top of the A set model state. The 
total ground state degeneracy for this situation has been 
studied in the literature: see Ref.[l^ for the Laughlin and 
Moore- Read states and Ref. [2^ for the Read-Rezayi and 
Ardonne-Schoutens series of non-abelian FQH states. 

For the Laughlin states the details are as follows. The 
TV-particle Laughlin state is realized on a total of N^f, + 1 
Landau orbitals, corresponding to = m{N ~ 1) flux 
quanta. The Laughlin state for ua particles would need 
— Tn{nA — 1) flux quanta; we thus have an excess 
flux of AiV^ = N^~Nf = m{N - ua). With the Laugh- 
lin gauge argument this corresponds to the presence of 
riqh = AiV^ quasi-holes over the groundstate. According 
to Ref. 13 each of the quasi-holes has a number of + 1 
effective orbitals to choose from, with bosonic counting 
rules (meaning that two or more quasi-holes can be in the 
same effective orbital). This gives a number of quasi-hole 
states equal to 

(ua + 1) + riqh ~ 1 
leading to the following upper bound to the entropy Sa 



S 



bound 



= In 



N^ + l^{m~l){nA-l) 



riA 



(4) 



We remark that this expression has a clear interpretation 
in terms of exclusion statistics: the counting factor in 
Eq. (ID) gives the number of ways ha particles can be 
placed in + 1 orbitals, in such a way that a particle 
placed in a given orbital I excludes particles from orbitals 
V with |; - l'\ < to. 

In a 1 /N expansion we find (assuming ua ^ N) 
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abound 



1 TO — 1 
N TO 

1 TO — 1 



2m?- 
-©(l/TV^) 



UA^UA - 1) 

nA{nA — 1)[2to -f (ua — 1)(to -f- — 4)] 



(5) 



The particle entropy reaches a maximum for ua = 
A^/2. For this case our Eq. ^ gives, in the limit of large 

S„^=N/2 < N[{m + 1) ln(TO + 1) - TO ln(TO)]/2 . (6) 
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This bound is sharper than a bound recently presented 
in Ref . [13, which gives a larger coefficient for the linear- 
in-A^ behavior. 

For the fermionic Moore-Read states a.t i/ = 1/m, with 
771 ~ 2,4,..., we can reason in a similar way, with now 
Nci, — m{N — 1) — 1. As for the Laughlin states we have 
an excess flux of AA^^ = — = 'm{N — nA) but 
now the number of quasi-holes is twice this number due 
to the fact that the fundamental quasi-holes correspond 
to half a flux quantum. Thus, riqh = 2AN^. We now take 
from Ref. Ul the following result for the total quasi-hole 
degeneracy 

^ f «qh/2 \finA- F)/2 + 77qh \ 

F=nA mod 2 ^ ^ ^ "'^^ ^ 

This gives us an upper bound 5^°^"'^ as before. 

Putting m — 2, one easily checks that 5^°'^"'^ coincides 
with S\ for HA — 2. In a expansion, the leading 
deviation from 5"^ is a 3-body term at order 

Si - 5^"' = -^l^AiuA - l){nA -2) + ... (8) 

This result nicely illustrates the fact that the leading cor- 
relations in the m = 2 Moore-Read state have a 3-body 
character: the wave-function vanishes if at least three 
particles come to the same position. 

For m ^ 2 the leading correlations do have a 2-body 
character, as for the Laughlin states. 

Si - ST-^ = L!!l_lnA{nA - 1) + . . . (9) 

Inspecting the particle entanglement at ua = N /2 and 
for N large, our bound implies that for the m — 2 Moore- 
Read state 

SnA=N/2 < 1.044iV . (10) 

This bound is reduced from the Fermi bound Sa ~ 
iV(41n4 - 31n3)/2 ~ 1.125 iV, but it is larger than the 
bound for the m = 2 (bosonic) Laughlin state, which 
has asymptotic form 7V(31n3 - 21n2)/2 - .955 7V. This 
indicates that, at equal filling v — 1/2, the particles in 
a Moore-Read state are more entangled than those in a 
Laughlin state. 

The quasi-hole counting rules for the order-A; clustered 
spin-polarized (Read-Rezayi) and spin-singlet (Ardonne- 
Schoutens) states are all known in the literature. They 
can be used to generalize the upper bounds on particle 
entanglement entropy given in this subsection to these 
more intricate non-abelian FQH states. 

C. Bosonic quantum Hall states 

We briefly comment on the case of bosonic FQH states. 
The realization that a rapidly rotating Bose gas may 



eventually enter a regime of bosonic quantum Hall states 
motivates the theoretical study of the effects of bosonic 
statistics. 

We consider bosonic Laughlin states at filling fraction 
1^ = ^ with 777 = 2, 4, . . .. The naive upper bound to the 
the entropy associated to placing ua bosons in + 1 
orbitals is 

^5 = In ( ) (11) 

The expression for 5*^°""^ remains unchanged, giving the 
following leading correction in a 1/N expansion 

^B_^bound^^„^(„^_l)^^^^ (12) 

For a bosonic Moore-Read state, with filling fraction 
v = 1/m with jn = 1,3,..., the leading 1/A^ correction 
becomes 

Sl-S\°^^'' = ^^^nA{nA^l) + ... (13) 

iV 777 

In the case m — I the leading correlations have 3-body 
character, leading to the vanishing of the leading 1/A^ 
correction. 

D. Numerical results 

In deriving the upper bound 5^°""'* we relied on the 
fact that a certain number of eigenvalues of the reduced 
density matrix vanish. The bounds would be exact if 
all non-zero eigenvalues were equal, but since they are 
not the bounds overestimate the actual values for the 
entropies. 

Fig. [1] plots the eigenvalues for the ua = 2-particle 
reduced density matrix for = 9 particles on a sphere 
in the m = 3 Laughlin state, for which the single par- 
ticle angular momentum is L = 12. The horizontal axis 
represents the degeneracy 2L2 + 1 of the eigenvalues, in 
descending order. The eigenvalue at L2 = 2L — 1 = 23, 
with degeneracy 47, vanishes; the non-zero eigenvalues 
show some scatter around an asymptotic value. Due to 
this scatter the entropy S — 5.509 is somewhat lower 
than the upper bound 5'^°"'^'^ = 5.533. 

An important difference between the 777, = 3 Laugh- 
lin and the 777, = 2 Moore-Read states is the absence of 
vanishing eigenvalues for the 2-particle reduced density 
matrix. The eigenvalue distribution shown in Fig. [2] il- 
lustrates this point. 

In the m = 2 Moore-Read state, there are vanishing 
eigenvalues in the reduced density matrix of ha > 3 par- 
ticles. The number of nonzero eigenvalues predicted by 
Eq. ([7]) agrees with numerical results. For example, for 
77^ = 3 and = 10 particles there are 770 nonvanishing 
eigenvalues, in agreement with Eq. ([7|). 

In Figs.[31l3]we compare numerically computed particle 
entanglement entropies with the bounds derived above. 
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E. Corrections to 5^4°""'^ due to eigenvalue spread 
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FIG. 1: (Color online) Eigenvalues for the 2-particle reduced 
density matrix, plotted against their multiplicities, for = 9 
particles in the m = 3 Laughlin state. 
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FIG. 2: (Color online) Eigenvalues for the 2-particle reduced 
density matrix, plotted against their multiplicities, for A'^ = 12 
particles in the m = 2 Moore-Read state. 
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FIG. 3: (Color online) Entanglement entropy for n,4 = 2 
and ua = 3 particles for the m = 3 Laughlin state. Dots are 
numerical exact values, the dotted line represents a.nd the 
solid curve is the bound S'a"""'^. 
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It is interesting to consider in some detail the devia- 
tion between the bounds 5'^°""'^ and the actual entropies 
computed numerically. As mentioned above, this devia- 
tion arises from the fact that the non-zero eigenvalues of 
the reduced density matrices are not all equal. 

To estimate the effect on Sua W] of the spread in the 
non-zero eigenvalues, we do a rough modeling of the 
eigenvalue distribution (Fig. [1]) of PnA=2 for the m = 3 
Laughlin state. For this case the number D of non-zero 
eigenvalues \s D ~ {3N — 4)(3A^ — 5)/2. If these nonzero 
eigenvalues were all equal (to 1/D), the entanglement en- 
tropy would have the maximum value Sd = InZ), which 
is the predicted upper bound ([4]). We now take into ac- 
count the deviations from l/D, guided by Fig. [TJ with 
the following toy distribution: we take Dq out of D of the 
eigenvalues to be equal to a/ D, with a > 1, while the rest 
of the eigenvalues are at value (3/D, (3 < I, such that the 
sum of eigenvalues is unity. Assuming Dq /D <C 1 leads 
to 



S ~ In D — {alna 



1 



(14) 



Guided by the eigenvalue distributions in Fig. [TJ we as- 
sume that Dq/D is of order 1/A^; for concreteness we 
put _Do equal to the multiplicity of the largest eigenvalue, 
which is 6A^ — 7. Taking a between 1.2 and 1.5 (as ob- 
served for the largest available Laughlin wavefunctions) 
gives a 1/A^ correction in the entropy with coefficient in 
the range 0.03 ~ 0.14. 

Fitting the difference S^ — SnA=2 to a form a/N+b/N"^ 
gives a coefficient a ~ 1.38. The vanishing eigenvalues 
account for a = 4/3, see Eq. ([5]), and we see that the 
remaining difference Sa ~ 0.05 is consistent with the 1/A'^ 
correction due to the spread in the non-zero eigenvalues. 

We made similar estimates for a the to = 2 Moore- 
Read state with up to TV = 12 particles, where the eigen- 
values are all non-zero and S']^"^"^ agrees with S^. In 
this case the deviation between data and bound show a 
1/N dependence with a coefficient of about 0.14. 

These considerations are of some general interest, as 
they make the point that a 1/A'^ expansion of particle 
entanglement entropies are indicative of correlations in 
a many-body state. In the concrete case studied here, 
the sizeable value of the leading 1/N correction in the 
Laughlin state indicates strong 2-body correlations, while 
the small value for the Moore-Read state indicates the 
absence of such correlations. 



Reduced density matrices and correlation 
functions 



FIG. 4: (Color online) Entanglement entropy for n,4 — 2 
and Ha = 3 particles for the m = 2 Moore-Read state. Dots 
are numerical exact values, the dotted line represents Sa and 
the soUd curve is the bound S*^""""^. 



Since the n^-particle reduced density matrices 
are obtained by integrating out all but ua of the par- 
ticles, one expects these matrices to be related to the 
n^i-particle correlation functions. In this subsection we 
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FIG. 5: (Color online) Two-particle correlation as a function 
of distance in units of the magnetic length, for the m = 3 
Laughlin state with A'' = 7, 8, 9 particles and for the m = 2 
Moore-Read state with N = 12. 

study this relation for the case ua = 2. In particular, the 
eigenvalue distributions of p„^=2 in Figs. [TJ [21 although 
discrete, are reminiscent of the well-known two-particle 
correlation functions (72 ('') for Laughlin and Moore-Read 
states. We will show that the eigenvalue distributions 
are in fact very closely related to the correlation func- 
tions ~ the eigenvalue distribution function is a kind of 
discretized version of 52 (r) ■ 

The two- particle correlation function g2ir) is conven- 
tionally defined as 



(15) 

where N is the number of particles and n is a density, 
which is chosen such that 52 (^) 1. 

We express the 2-particle reduced density matrix p2 
on a sphere in a basis of polar spherical coordinates. Be- 
cause of the rotational symmetry it should be a function 
of the angular distance 9 between the two particles. In 
Appendix [X] we show that p2 (0) can be written in the 
form 

2/ I 1 

P2{0) = Y.Xi^Me). (16) 
I 

In this expression. A; is the eigenvalue with multiplicity 
21 -\- 1, corresponding to the total angular momentum 
of the two particles equal to /. The functions Ri{6) are 
explicitly given in Eq. (jA12p . 

Since the distance r between two particles is simply 
equal to r = R6, with R the radius of the sphere, the 2- 
particle correlation function g2 (r) is directly proportional 
to p2 {0) . Through Eq. ([16]) it is expressed as a transfrom 
from I space to 6 space, with basis functions Ri{9). In 
Fig. [HI we show some curves for P2{9), they agree with 
known results. 



As illustrated in the inset to Fig.[5l the basis functions 
Ri{9) have a peak structure, with the position of the peak 
depending on the total angular momentum I. Large val- 
ues of / correspond to small angular distances and vice 
versa. This is easy to understand from the following clas- 
sical picture. When two particles with angular momenta 
L have total angular momentum 2L, the corresponding 
vectors Li and L2 should point a the same direction. 
The angular momentum vector points at the position of 
the particle on a sphere, therefore two particle should be 
close to each other. On the other hand, if total angular 
momentum is zero then the angular momentum vectors 
should point into opposite directions. This means that 
particles are placed at the opposite sides of the sphere. 

The fact that the Ri{9) are localized functions, peaked 
at 9 values monotonically decreasing with indicates 
that g2 {r) curve is simply a continuous form of the A; ver- 
sus descending- (2Z + 1) curves of Figs. [T] and [5) The sim- 
ilarity between the discrete A; and the continuous 32 (j") 
curves is not accidental. 

At small distances Ri{9) oc (0^)^^"'. For a Laughhn 
state the lowest value of 2L — Z is m, thus P2{9) oc 9'^'"\ 
This behavior is a direct consequence of the vanishing of 
eigenvalues with the largest multiplicities. For the m = 2 
Moore- Read state the lowest value of 2L — / is 1 because 
there are no vanishing eigenvalues. Therefore at small 
distances P2{9) oc 9^. 

Of course, our observations on the 2-particle correla- 
tions agree with known results; our main point has been 
to stress the intimate relation with the eigenvalue distri- 
bution of the 2-particle reduced density matrices. 



IV. SPATIAL ENTANGLEMENT AND 
TOPOLOGICAL ENTROPY 

We now turn to dividing the Landau level orbitals into 
two blocks and calculating the entropy of entanglement 
between them. In a previous publication^ we used this 
scheme to extract the topological entropy of the Laughlin 
state. Here, we mainly focus on the (m — 2) Moore-Read 
state. After reviewing the topological entropy 7 and the 
total quantum dimension V, especially in the context of 
the Moore- Read state (jIV Ap . we detail some issues with 
taking the thermodynamic limit necessary for extract- 
ing 7 from numerical data pVBp . and then present our 
numerical results pV Cp . We also present observations 
on the spectral structure of the reduced density matrices 

(irvDi) . 



A. Topological entropy for the Moore-Read state 

For spatial partitioning of many-particle states, the 
general rule ("area law") is that the entanglement en- 
tropy scales as the size of the boundary between the 
A and B blocks. Subtle information about the nature 



7 



of the many-particle state can be provided by the pres- 
ence or absence of logarithmic corrections, values of co- 
efficients, or subleading terms in this basic relationship. 
For topologically ordered states in two dimensions, the 
following theorem has been presented recently^^'^^ con- 
cerning the scaling of entanglement entropy between spa- 
tial partitions. If L is the length of the boundary be- 
tween the two blocks, the entanglement entropy scales as 
Sa — ctL — 7 -f 0(L~^). As usual the scaling law apphes 
to situations where A is large and the total system is 
infinite. The subleading term 7 is called the topological 
entanglement entropy. The striking result of Refs. Illlfl^ 
has been that it can be expressed as the logarithm of a 
quantity 2? known as the total quantum dimension of the 
topological field theory describing the topological order 
of the state. The total quantum dimension is given by 



V 



(17) 



where the d^'s are the quantum dimensions of the in- 
dividual sectors making up the topological field theory. 
These quantum dimensions are set by fusion rules of the 
fundamental anyons in the field theory. 

The topological field theory for a = 1/m Laughlin 
state has a fundamental anyon (of fractional charge q = 
e/m), which generates m abelian sectors. The quantum 
dimension di is unity in all sectors so that V = ^Jm, 
for the V = 1/m Laughlin state. For m = 3 this gives 
7 = lnV3 ~ 0.55. 

For states with non-abelian quasiparticles, the situ- 
ation is more interesting because some anyon sectors 
contribute c?i > 1. Details for some examples have been 
provided in Refs. 

In particular, for the m — 2 
Moore-Read state, there are six sectors (two each of 
quasiparticles denoted by /, a, ip) which contribute 
di = \, da = V^, = 1' leading io V — \/8 and 
7 ~ 1.04. The non-abelian nature shows up in the fact 
that 7 is larger than In-v^G, six being the degeneracy of 
the m — 2 Moore-Read state on the torus. 

A compact general expression for the total quantum 
dimension for a Read-Rezayi state with order-/c clustering 
and at filling fraction v = k/{kM + 2) is (see also Ref. 
24 



VnR[k,M] 



V(fc + 2)(fcM + 2) 
2sin(7r/(fc + 2)) ' ^ ^ 

It includes the Laughlin states (fc = 1, M = m — 2) 
and the Moore-Read states {k = 2, ill = m — 1) as 
special cases. For a general spin-singict non-abelian 
FQH state with order k clustering and filling fraction 
1/ = 2k/{2kM + 3), the result is 



(fc + 3)V(2fcAi" + 3) 



(19) 



givmg 7 
M = 1) at 



16cos(7r/(fc + 3))sin3(7r/(fc + 3)) ' 
1.62 for the paired spin-singlet state {k — 2 



In general, the M-dependence of these expressions for 
total quantum dimensions is linked to the ground state 
degeneracy in torus geometry. Denoting the latter by 
=i^[fc, M] we have the relation 



V[k,M] ^V[k,Q] 



#[fc,o] 



(20) 



4/7. 



The conformal field theories underlying the states at 
M = are of Wess-Zumino-Witten type {SU{2)k for the 
RR states and SU{3)k for the AS series) and the quan- 
tities 0] can be expressed in the modular 5-matrix 
for these WZW models. 



B. 7 from sphere calculations: extrapolations 

First of all, we note that, since our degrees of 
freedom are ordered essentially one-dimensionally, we 
cannot use one of the two-dimensional schemes pro- 
posed prcviousl y^^'^^1^^ in which an appropriate addi- 
tion/subtraction of the entanglement entropies of sev- 
eral regions cancels the boundary parts of the entropy 
(Sa olL — 7) leaving the subleading term 7. With the 
orbital degrees of freedom on a sphere, we can choose 
only regions corresponding to disks, concentric rings, and 
combinations thereof. Any combination of entropies of 
disk- and ring-like regions that cancels out the boundary 
terms also unfortunately cancels out the 7 term. 

We are thus led to using directly the scaling law, 
Sa '"^°°> aL — 7. Our choice of block A as the first I a 
orbitals, extending spatially from one pole out to some 
latitude, corresponds to a disk-shaped block only in the 
thermodynamic limit. The block area is proportional to 
the square of ^/Ta while its boundary is proportional to 
y^lA{N,j, + 1 — Ia)', these are equivalent only in the same 
— > 00 limit. One way to numerically access the ther- 
modynamic limit is to take the entanglement entropy of 
I A orbitals with the rest, for accessible wavefunctions of 
various sizes N, and then take the N ^ 00 limit. The 
Si^{N 00) versus ^/Ia points thus obtained should 
then follow a linear curve at large Ia^ whose vertical in- 
tercept gives the topological entropy. Results following 
this procedure were provided for the v = 1/3 Laugh- 
lin state in our earlier paper'^ here we will focus on the 
Moore-Read state. The extrapolation of S*;^ (A^) values to 
the thermodynamic limit is a tricky issue. We therefore 
discuss the extrapolation in some detail here, providing 
some general results. 

We are interested in the function Si^ (x) , where x — 
1/N. We have access to Si^(xi) at several integer val- 
ues of A^, and would like to estimate Si^(0). For each 
dataset that we have access to (each Ia; both Laugh- 
lin and Moore-Read), we note the following: the Si^{xi) 
versus Xi values form a monotonic curve and this curve 
gets flatter (slope magnitude decreases) with decreasing 
X. Two examples can be seen in the inset to Fig. [51 
In other words, the first and second derivatives of the 
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UJ = 2.0 


3.46601 


5.21183 


4.69023 


4.77647 


3.80408 


4.92355 


4.74756 




4.01821 


4.84476 






4.16743 








oj = 3.0002 


3.46601 


4.56058 


4.51976 


4.55722 


3.80408 


4.53643 


4.59865 




4.01821 


4.56761 






4.16743 








OJ = 2.5926 


3.46601 


4.74638 


4.58027 


4.65108 


3.80408 


4.65109 


4.65108 




4.01821 


4.65108 






4.16743 











TABLE I: Extrapolation using the BST algorithm, using 
three different u values for the same initial dataset (first col- 
umn), corresponding to Sij^=]_i{N) values for N = 12,14,16, 
18. In the last case, uj has been tuned to ensure that the 
sets obtained after first and second iterations converge to the 
same value, i.e., the lowest elements of the second and third 
columns are the same. 

5;^ {x) function have the same sign and neither deriva- 
tive changes sign. 

Motivated by the above observations, we provide the 
following result. Assuming only that the signs of the first 
two derivatives of the Si^ {x) function are the same and 
that the signs remain unchanged until a; = 0, we have: 

1. The value Sq = Sij^ (xg) corresponding to the small- 
est value a;o of the available strict lower 
(upper) bound for S*;^ (0) if the {x) is negative 
(positive). 

2. The intercept found by connecting the 5;^ (a;) cor- 
responding to the smallest two Xi values (xq, xi), 
namely 

5i Si^{xf)) { 1 — — ) + Si^{xi) { — — — ) , 

V xi-xaj \xi-xoJ 

is a strict upper (lower) bound if S'l^ (x) is negative 
(positive) and S['^{x) is positive (negative). 

The limits Sq and Si thus obtained give us conserva- 
tive bounds for the required entanglement entropies in 
the thermodynamic limit, Sij^{N —f oo). To obtain a 
sharper extrapolation, one can use various polynomial 
extrapolations and take the average, as done in our ear- 
lier work.^'' Here, we improve the extrapolation by us- 
ing the extrapolation algorithm of Bulirsh and Stoer 
(BST algorithm), based on rational polynomial fraction 
approximations 

The BST procedure involves successive transforma- 
tions of the original dataset, leading to transformed 



datasets successively smaller by one element. If the free 
parameter (w) of the algorithm is chosen properly, the 
successive sets will be more and more convergent and will 
eventually converge to the x = value. The parameter 
LD is chosen to optimize convergence; how precisely this is 
done is a non-trivial issue, especially for large datasets. 

In the present case where the initial dataset consists of 
bipartite entanglement entropies, it is particularly desir- 
able to be able to use a small number of initial Si^{N) 
values, corresponding to the largest available N . This 
is because the symmetry S*;^ ~ Sn^+i~Ia makes the 
Si A (N) data meaningless for N oo extrapolation when 
I A gets close to Nci,/2 — \{mN — S). In other words, the 
Sij^[N) values used should not be too close to the peak 
at the midpoint of the respective Si^ (N) versus I a curve 
(e.g., Figure 1 of Ref. 10). This restricts us to the largest 
few N values. 

In Table U we demonstrate some possible choices of uj 
starting from four S';^(A^) values. The top example uses 
Lu = 2, which is a common choice used for extrapolation- 
based numerical integration algorithms. In the second 
example, uj is tuned to give the "most converged" set 
after the first iteration, by choosing uj to minimize the 
standard deviation of the second column. We find that 
this procedure does not always lead to great stability for 
the subsequent iteration results. 

In the last example, we have tuned uj to give the same 
last value for the sequences obtained after the first and 
the second iterations. This also guarantees the third (fi- 
nal) iteration to give the same value (Table |T|, which 
indicates that the convergence is very good. This proce- 
dure has the added advantage that actually only the last 
three points of the initial dataset are used to determine 
the values tuned to be equal. We thus get a nicely con- 
verged estimate based only on three initial points, which 
is highly desirable as explained above. This is therefore 
our method of choice for tuning the free parameter uj in 
applying the BST algorithm for our problem. 

As a side note, we remark that if the largest available 
entanglement entropies are not known exactly (e.g., if 
they are calculated from approximate wavefunctions, or 
using the desnity matrix renormalization group (DMRG) 
technique), it might be more prudent to use a prescrip- 
tion for choosing uj that actually uses more than the min- 
imal (three) numbers that we have used here. In the 
cases reported here and in our earlier workji^ we have 
entanglement entropies calculated from numerically ex- 
act wavefunctions, so we have no precision issues to worry 
about when choosing the number of Si ^ (N) values to use 
for the N ^ oo extrapolation. 

Once the extrapolation value is (uniquely) determined 
by the algorithm outlined above, we need an estimate of 
the uncertainty. We can use the bounds (5*0, S'l) derived 
earlier to obtain a conservative error estimate: 

min(|5i - 5bst|, I'5'bst - ^ol) , 

where S'bst is the extrapolated value obtained by the 
BST algorithm. 




Sq. root of # orbitals in block, / 



FIG. 6: (Color online.) Entanglement entropies in Moore- 
Read state wavefunctions, extrapolated to the thermody- 
namic limit. Dashed line is a fit to —7 -|- CivTa, with some 
points dropped. Inset plots 5;^ against 1/A'' for various fixed 
Ia. 
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FIG. 7: (Color online.) Density matrix eigenvalues in de- 
creasing order for orbital partitioning with Ia orbitals in a 
block for = 9, m = 3 Laughlin state. 



also some error propagated from the extrapolation uncer- 
tainty. Using the conservative uncertainty estimate pro- 
posed in the previous subsection gives us a more conser- 
vative and more rigorous error estimate, 7 ~ 0.51 ± 0.25. 



C. Numerical results 

Moore-Read state. Fig. [5] shows results of numerical 
calculations for the v = 1/2 Moore- Read state. We 
used exact wavefunctions up to = 18 particles. These 
wavefunctions were obtained by diagonalizing LP' in an 
Lz —Q Hilbert space spanned by the "squeezed states"—. 
After numerically obtaining the entanglement entropies 
Sij^ (TV) from these wavefunctions, we obtain estimates 
and uncertainties for the N ^ 00 extrapolations by the 
procedure outlined in the previous subsection. The re- 
sulting data are plotted in Fig. [51 

The linear S*;^ versus \/Ta behavior is expected only 
for large I a', however our lavge-lA points have the greatest 
uncertainty. For estimating the topological entropy, we 
therefore make linear fits after discarding to 5 of the 
smallest-?A points and/or to 2 of the largest-/^ points. 
This results in estimates of 7 (magnitude of the vertical 
intercept) scattered between 0.85 and 1.35. The error 
propagated into our 7 estimate from our extrapolation 
uncertainties is ~ 0.3, larger than that obtained from this 
scatter. With all this we arrive at the result 7 ~ 1.1±0.3, 
quite consistent with the expected value of 7 ~ 1.04. 

Laughlin state. We used the well-defined procedure 
of the previous subsection to revisit our previous esti- 
mate of the topological entropy for the = 1/3 Laugh- 
lin statejiS To get the extrapolated Si^ , we now use the 
BST estimates rather than doing several polynomial fits. 
Dropping to 4 of the smallest-?^ points and to 2 
of the Irgest-^A points leads to 7~0.51±0.14, consistent 
with the previously reported estimate (0.60 ± 0.15) and 
with the expected value 7 « 0.55. The error estimate 
reported in Ref. [lOl only took into account this varia- 
tion, due to dropping various number of points. There is 



D. Eigenvalue distribution for reduced density 
matrix 

In Fig. [71 we show the largest eigenvalues of reduced 
density matrices obtained by orbital or spatial partition- 
ing. The eigenvalues are ordered according to decreasing 
magnitude and plotted on a log scale; the resulting curves 
are roughly linear, suggesting a roughly exponential de- 
cay of the eigenvalue distribution function. 

It is interesting to note the complete dissimilarity of 
this eigenvalue spectrum compared to the particle parti- 
tioning case discussed earlier, e.g., Fig. [H It would also 
be interesting to put our observations in the context of 
the spectra of reduced density matrices of many-body 
systems in general. Reduced density matrices for spa- 
tially connected blocks have been studied previously in 
the context of the convergence of the DMRG algorithm; 
an overview is available in section III-B of SchoUwock's 
DMRG review From our numeric data, It is difficult 
to say whether or not the decay of the eigenvalue distri- 
butions is slower than exponential. 



V. CONCLUDING REMARKS 

We have presented a detailed study of the entan- 
glement entropy in abelian and non-abelian quantum 
Hall states, taking a paradigmatic example of each, the 
ly = 1/3 Laughlin state and the = 1/2 Moore- Read 
state. 

For entanglement between subsets of particles, we have 
demonstrated the effects of particle-particle correlations 
in the deviation of the entanglement entropies from an 
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upper bound set by fermionic statistics only. We pre- 
sented a close upper bounds for particle entanglement 
entropies for both bosonic and fermionic states. For a 
particular case ha — N/2, i.e. entanglement of the one 
half of the system with the other, the presented upper 
bound is notably lower than a similar estimate in Ref. 
We showed that the distribution of the eigenvalue spec- 
trum of the two-particle density matrix can be directly 
related to a two-particle correlation function 52 (f) ■ This 
suggests that the eigenvalue distribution of n^-particle 
reduced density matrices can be related to n^-particle 
correlation functions. 

For entanglement between spatial regions, indexed by 
orbitals on a spherical geometry, we have provided a thor- 
ough discussion of our efforts to optimize the procedure 
for extracting the so-called "topological entropy" from 
finite-size data. Using this procedure, we calculated the 
topological entropy for the Moore-Read state, and found 
our numerical result to be consistent with the expected 
value In VS. 

Our results on Laughlin and Moore-Read states indi- 
cate that the computation of topological entanglement 
entropies can be used for diagnosing topological order in 
FQH states that are only available for a limited number of 
particles. Clearly, the practical use of this method hinges 
on the accuracy by which a value for 7 can be extracted. 
In this context, it could be particularly interesting to use 
entropies 5/^ [N] that are generated in DMRG studies of 
FQH states. 

It is interesting to note the differences between the 
particle and spatial partitioning cases. In the first case, 
correlation effects decrease the particle entanglement en- 
tropy Sua ■ In the latter case, correlation effects increase 
the spatial entanglement entropy; for example, a Wigner 
crystal would have very high entanglement entropy be- 
tween two spatial partitions. These differences are partic- 
ularly dramatically manifested in the difference between 
the spectral structure of the reduced density matrices for 
the two cases. In the particle entanglement case, the 
eigenvalues tend to be of the same order, with relatively 
small fluctuations around some average value. In the 
spatial case, the eigenvalue distribution function has an 
exponential-like form, and only a small fraction of the 
eigenvalues contribute to the entanglement entropy. 
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APPENDIX A: 2-PARTICLE REDUCED 
DENSITY MATRIX AND CORRELATION 
FUNCTIONS 



In this appendix we derive Eq. (|16p . The two-particle 
reduced density matrix may be written as 



P2 = y^^XiTn\l,rn){l,', 



(Al) 



where \l,rn) is a state of two particles with total angular 
momentum I and projection m. The rotational symmetry 
of the sphere (p2 must be invariant under rotations) indi- 
cates that the eigenvalues Xim do not depend on the pro- 
jection m, i.e, are grouped in SU{2) multiplets; A;„i = A/. 

To get the correlation function, we need to compute the 
density matrix in a basis of polar spherical coordinates. 



P2(6'i,6'2,01,02j 



71,1^2,01,021^2 |o'l,6'2,01,02, 



(A2) 

where 0i and 0^ are coordinates of particles on a sphere. 
Rotational symmetry requires that p2 depends only on 
the angular distance 9 between particles: 



cos = cos( 01 — 02 ) sm Hi sm 6'2 + cos Oi cos 6^2 



(A3) 



Now the unnormalized wavefunction of two particles 
with total angular momentum I can be written as^ 



= Y[{a*u,+f3*Viy{uiV2~U2Vi) 



2L-1 



(A4) 



1=1,2 



where L is the angular momentum of each particle, Ui = 
cos ^6**^'^^, Vi = sin ^e~*"*'/^, and a, (3 fix the center of 
mass of the two particles. The same wavefunction can 
also be written as 



,0i,02|i?(a,/3)|;,/) 



(A5) 



where D{a,(3) is a rotation operator. Let us write jfi) = 
D{a,f3) \l,l), which is a familiar coherent state. Using 
the identity 



21 + 1 
in 



J dn\fi){fi\ = ^|;,m)(/,i 



(A6) 



one obtains 

^1(01,02,01,02 |;,m)| 

m 

21+1 



rff! 1(01,02,01,02 (A7) 



The normalization of (0i, 02, 0i, 02 |^^) is equal to2^ 
167r2(n)2 {2L - 1)\{2L + I + l)\ 



N, 



(21 + 1) 



(2L + 1)P 



(A8) 
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We need to compute 



we finally write the 2-body correlator as 



ii = J dn 1(01,02,01,02 

^ A in 



1=1,2 



(A9) 



Defining a = \u\V2 — 1*2^11^ = (1 — cos0)/2 and h 
\u\ui + v-2v\ P = (1 + cos 0)/2, one finds 



{21 



47r(;!)2 



(21 + 1) 

One can see that Ii is normalized as 



^a^"^ 2Fi{-l, 1 + 1,1, ~b/a) . (AlO) 



d9sin0Ii{e) 



27r 



(All) 



Introducing the function 



Ri{e) 



Ni 



(2^+1)!' 



At:{2L~1)\{2L + 1 + 1)\ 



a'"^ 2Fi{-l, 1 + 1,1, -b/a) 

(A12) 



(A13) 



APPENDIX B: ENTANGLEMENT ENTROPIES; 
NUMERICAL RESULTS WITH ORBITAL 
PARTITIONING 



In tables |TT] and IIIII we list orbital-partitioning entan- 
glement entropies calculated using the numerical wave- 
functions. For each wavefunction (each colmun), the en- 
tanglement entropies are only listed up to their maximum 
value, typeset in bold, because the values after this are 
determined by the symmetry Si^ = S^^+i-i^. The ex- 
trapolation procedure of Sec. IIV Bl involves extrapolating 
each row of numbers to iV ^ cx3. 

Note that, in each column (for a particular TV), the 
omitted part after the midpoint (in bold) is a decreasing 
function of I a and thus does not give useful information 
about the thermodynamic limit of Si^. In the extrapola- 
tion, it is therefore important to avoid values from these 
parts of the table. We therefore restrict ourselves to val- 
ues which in tables [TTl and [1111 are to the right of (i.e., 
above) the diagonal line through the midpoint numbers 
typeset in bold. 
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Ia 


N = 4 


N = 6 


iV = 8 


N = 10 


N =12 


Af = 14 


N = 16 


N =18 


1 


0.636514 


0.673012 


0.682908 


0.686962 


0.689009 


0.690186 


0.690923 


0.691416 


2 


1.09861 


1.15777 


1.23519 


1.24961 


1.25796 


1.2663 


1.27095 


1.27385 


3 


1.09861 


1.49971 


1.65513 


1.7002 


1.72587 


1.74591 


1.75859 


1.76716 


4 




1.76712 


2.0355 


2.11077 


2.15836 


2.19257 


2.21444 


2.22959 


5 




1.88152 


2.31807 


2.44895 


2.52664 


2.58123 


2.61639 


2.64089 


6 






2.48295 


2.70298 


2.82429 


2.90754 


2.96113 


2.99855 


7 






2.54282 


2.88348 


3.06178 


3.18053 


3.257 


3.31045 


8 








2.98996 


3.24156 


3.40254 


3.50653 


3.57927 


9 








3.02558 


3.36728 


3.57792 


3.71402 


3.80922 


10 










3.44165 


3.71089 


3.88358 


4.00428 


11 










3.46601 


3.80408 


4.01821 


4.16743 


12 












3.85949 


4.12054 


4.30135 


13 












3.87774 


4.19219 


4.408 



TABLE II: Orbital-partitioning entanglement entropies <S'(^[A'^] for the i/ = 1/2 Moore- Read state. 



Ia 


N = 3 


N = 4 


N = 5 


N = 6 


N = 7 


N = 8 


N = 9 


iV = 10 


1 


0.682908 


0.673012 


0.666278 


0.661563 


0.65811 


0.655482 


0.653418 


0.651757 


2 


1.00424 


1.05492 


1.07339 


1.0822 


1.08707 


1.09006 


1.09202 


1.09337 


3 


1.27703 


1.3944 


1.44547 


1.46998 


1.48491 


1.49466 


1.50147 


1.50648 


4 


1.27703 


1.59387 


1.70596 


1.76262 


1.79672 


1.81933 


1.8353 


1.84717 


5 




1.66184 


1.87354 


1.97796 


2.03976 


2.08068 


2.10962 


2.13113 


6 






1.9523 


2.1202 


2.21872 


2.28347 


2.32919 


2.36314 


7 






1.9523 


2.20091 


2.3448 


2.43838 


2.50425 


2.55307 


8 








2.22768 


2.42555 


2.55286 


2.64196 


2.70783 


9 










2.46451 


2.63139 


2.74713 


2.83233 


10 










2.46451 


2.67761 


2.82363 


2.93049 


11 












2.69293 


2.87356 


3.00467 


12 














2.89817 


3.05653 



TABLE III: Orbital-partitioning entanglement entropies Sij^[N] for the z/ = 1/3 Laughlin state. 
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